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Abstract 

Miniature Inverted-repeat transposable elennents (MITEs) are nunnerlcally predominant transposable elements in the rice 
genome, and their activities have influenced the evolution of genes. Very little is known about how MITEs can rapidly 
amplify to thousands in the genome. The rice MITE mPing is quiescent in most cultivars under natural growth conditions, 
although it is activated by various stresses, such as tissue culture, gamma-ray irradiation, and high hydrostatic pressure. 
Exceptionally in the temperate japonica rice strain EG4 (cultivar Gimbozu), mPing has reached over 1000 copies in the 
genome, and is amplifying owing to its active transposition even under natural growth conditions. Being the only active 
MITE, mPing in EG4 is an appropriate material to study how MITEs amplify in the genome. Here, we provide important 
findings regarding the transposition and amplification of mPing in EG4. Transposon display of mPing using various tissues of 
a single EG4 plant revealed that most de novo mPing insertions arise in embryogenesis during the period from 3 to 5 days 
after pollination (DAP), and a large majority of these insertions are transmissible to the next generation. Locus-specific PGR 
showed that mPIng excisions and insertions arose at the same time (3 to 5 DAP). Moreover, expression analysis and in situ 
hybridization analysis revealed that Ping, an autonomous partner for mPing, was markedly up-regulated in the 3 DAP 
embryo of EG4, whereas such up-regulation of Ping was not observed in the mP/ng-inactive cultivar Nipponbare. These 
results demonstrate that the early embryogenesis-specific expression of Ping is responsible for the successful amplification 
of mPing in EG4. This study helps not only to elucidate the whole mechanism of mPing amplification but also to further 
understand the contribution of MITEs to genome evolution. 
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Introduction 

Transposable elements (TEs) are DNA sequences that are 
capable of jumping from one genomic locus to another and make 
up a large fraction of eukaryotic genomes. More than 80% of the 
maize {^ea mays) and barley [Hordeum vulgare) genomes are 
composed of TEs [1], [2], and they constitute 35% and 14% of 
the genomes of rice [Oryza sativa) and Arabidopsis [Arabidopsis 
thaliand], respectively [3] , [4] . TEs are harmful to the host because 
their mobilities perturb genome stability, whereas they play greatly 
generative roles in genome evolution such as alternation of gene 
structure, change of expression pattern, and rearrangement of 
chromosome structure [5], [6]. 

TEs are classified into two groups according to their transpo- 
sition mechanisms: class I elements (retrotransposons) that 
transpose through a copy-and-paste mechanism via an RNA 
intermediate, and class II elements (transposons) that transpose 
through a cut-and-paste mechanism without undergoing an RNA 
intermediate. Class I elements easily attain tens of thousands of 
copies, whereas the majority of class II elements cannot amplify 
themselves to 50 copies at most. Unlike other class II elements. 



miniature inverted-repeat transposable elements (MITEs) have die 
capacity to amplify themselves to high copy numbers (hundreds or 
thousands) [7]-[9]. In the rice genome, MITEs are numerically 
predominant TEs [10], constituting 8.6% of the genome [11]. 
Because MITEs are too short (<600 bp) to encode any protein, 
their transpositions must depend on the proteins encoded by the 
autonomous elements. Well-studied MITEs are classified into the 
Stowaway and Tourist families, which belong to the Tel /mariner and 
PIF/ Harbinger superfamUies, respectively. Because MITEs are 
mainly deployed in gene-rich regions [10], [12] and affect 
adjacent gene expression [13], they are considered to play an 
important role in genome evolution. However, little is known 
about how MITEs attain high copy numbers. 

Miniature Ping (mPing) is the first active MITE identified in the 
rice genome [14]-[16]. Although MITEs are deployed in the 
genome at a high copy number, the copy number of mPing 
exceptionally remains at a low level in most rice cultivars: indiea 
and tropical japonica cultivars have fewer than 10 copies, and 
temperate japonica cultivars including Nipponbare have approxi- 
mately 50 copies [14]. The transposition oi mPing is suppressed in 
most rice cultivars, but, like other TEs, mPing is activated by 
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Author Summary 

Transposable elements are major components of eukary- 
otic genomes, comprising a large portion of the genome in 
some species. Miniature inverted-repeat transposable 
elements (MITEs), which belong to the class II DNA 
transposable elements, are abundant in gene-rich regions, 
and their copy numbers are very high; therefore, they have 
been considered to contribute to genome evolution. 
Because MITEs are short and have no coding capacity, 
they cannot transpose their positions without the aid of 
transposase, provided in trans by their autonomous 
element(s). It has been unknown how MITEs amplify 
themselves to high copy numbers in the genome. Our 
results demonstrate that the rice active MITE mPing is 
mobilized in the embryo by the developmental stage- 
specific up-regulation of an autonomous element. Ping, 
and thereby successfully amplifies itself to a high copy 
number in the genome. The short-term expression of Ping 
is thought to be a strategy of the mPing family for 
amplifying mPing by escaping the silencing mechanism of 
the host genome. 

exposure to various stress conditions such as gamma-ray 
irradiation [16], hydrostatic pressurization [17], cell culture [14], 
anther culture [15], and inhibition of topoisomerase II [18]. 
Introgression of distantly related genomes also causes mPing 
transposition [19], [20]. However, mPing is actively transposing 
without such stresses in the temperate japonica rice strain EG4 
(cultivar Gimbozu) under natural growth conditions, and its copy 
number is approximately 1000 copies [21]. This indicates that 
mPing has overcome the silencing mechanism or established a 
novel strategy for its amplification in the EG4 genome. In this 
sense, mPing in EG4 is an appropriate material to study the 
amplification of MITEs in plant genomes. 

The autonomous element Ping and its distantly related element 
Pong, which both belong to the PIF/ Harbinger superfamUy, provide 
two proteins required for mPing transposition. Both Ping and Pong 
have two open reading frames (ORFs), ORFl and ORF2 [22], 
[23]. The former encodes a Myb-hke DNA-binding protein, and 
the latter encodes a transposase lacking DNA binding domain. 
Transposase of most class II elements contains a conserved 
catalytic domain (DDE motif) and a DNA-binding domain [23], 
[24], whereas these domains are encoded separately by two ORFs 
in both Ping and Pong [22], [23]. The study of other members of 
the PIF/ Harbinger superfamily suggested that the Myb-like DNA- 
binding protein directly binds to the subterminal regions of the 
transposon in order to recruit the transposase [25]. Both Myb-like 
protein and transposase of either Ping or Pong or both elements are 
necessary for mPm^ transposition [22], [23]. 

In this study, we demonstrate that mPing is actively transposing 
in the embryo of EG4 during the period from the regionalization 
of shoot apical meristem (SAM) and radicle to the formation of the 
frrst leaf primordium (3 to 5 days after pollination, DAP) with the 
aid of developmental stage-specific expression of Ping. Our results 
provide important evidence for the amplification mechanism not 
only of mPing but also of other MITEs. 

Results 

Transpositions of mPing during gametogenesis 

Plants have acquired the silencing mechanism of TEs in germ 
cells. In Arabidopsis, for example, TEs are activated specifically in 
the vegetative nucleus of the pollen, and siRNAs from the 



activated TEs accumulate in the sperm cells [26] . On the basis of 
these results, Slotkin and colleagues proposed that siRNAs derived 
from TEs activated in the vegetative nucleus silence TEs in the 
sperm cells [26]. We conceived that mPing might overcome such a 
silencing mechanism in EG4. To confirm this hypothesis, we 
developed two Fi populations from reciprocal crosses between the 
mPffzg-active strain EG4 and the mPm^-inactive cultivar Nippon- 
bare, and investigated the transposition activity of mPing by 
transposon display (TD) analysis. Success of reciprocal crosses was 
confirmed by PGR analysis using locus-specific primers (Figure 
SI A). One of the results of TD analysis using two selective bases is 
shown in Figure lA; all 16 possible primer combinations were 
analyzed. The banding patterns of Fi plants were almost the same 
as those of EG4. The bands that appeared in all F] plants but not 
in the parental EG4 plant were derived from another parental 
Nipponbare plant (Figure SIB). Furthermore, the bands that 
appeared in only one of eight Fi plants but not in the parental 
EG4 plant are herein referred as de novo insertions. These bands 
were confirmed not to be PGR artifacts by sequence and locus- 
specific PGR analysis (Table SI and Figure S2). We detected 15.5 
de mm insertions per plant in the selfed progenies of EG4, whereas 
Nipponbare yielded no de novo insertions in the selfed progenies 
(Figure IB). This confirmed that mPing is active in EG4 under 
natural growth conditions but inactive in Nipponbare. If mPingwSLS 
specifically activated in the poUen of EG4, we could obtain de novo 
insertions only in the Fi plants from the Nipponbare/EG4 cross. 
However, we obtained de novo insertions in both Nipponbare/EG4 
and EG4/Nipponbare populations (Figure IB). Moreover, there 
was no significant difference in the number of de novo insertions per 
plant between the two Fi populations. This indicates that the 
activating factor(s) for the mPing transposition is present in both 
male and female gametes of EG4. 

Transpositions of mPing during ontogeny of EG4 plants 

We performed TD analysis of mPing using genomic DNA 
samples extracted from endosperm, radicle, and leaf blades of 
eight progenies (Si) derived from a single parental EG4 plant (Sq), 
and investigated the mPing transposition during ontogeny of rice 
plants (Figure 2A). One of the results of TD analysis using two 
selective bases is shown in Figure S3; all 16 possible primer 
combinations were analyzed. We examined de novo insertions in the 
same way as described above. Gonsequentiy, a total of 228 de novo 
insertions were detected. These insertions were divided into five 
groups (Figure 2B): (1) endosperm-specific insertions that appeared 
only in the endosperm sample, (2) radicle-specific insertions that 
appeared only in the radicle sample, (3) leaf-specific insertions that 
appeared only in one sample from the 1st to 3rd leaf blades, but 
not in the 4th and 5th leaf blades, (4) shoot-specific insertions that 
appeared in at least one sample of 1st, 2nd, and 3rd leaf blades, 
and in at least one sample of 4th and 5th leaf blades, and (5) 
radicle/ shoot-specific insertions that appeared in both radicle and 
leaf blade samples. These de novo insertions were confirmed by 
sequence and locus-specific PGR analysis (Table S2 and Figure 
S4). Numbers of each insertion obtained in this study are 
summarized in Figure 2G. Plant development is divided roughly 
into three successive phases: embryogenesis, vegetative phase, and 
reproductive phase. If m-P»z^ transposed in the SAM of the Sq plant 
during vegetative and/or reproductive phases, the de novo 
insertions would segregate according to MendeFs law among the 
Si progenies. We obtained no band that appeared in at least two 
S 1 progenies and was not seen in the So plant. This indicates that 
the transmissible insertion of mPing to the next generation seldom 
(or never) arises during the vegetative and reproductive phases. 
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Figure 1 . Transposition of mPing in reciprocal crosses between 
EG4 and Nipponbare. (A) Transposon display (TD) for mPing of the F, 
population from reciprocal crosses between EG4 and Nipponbare. One 
of the results of TD analysis using two selective bases is shown. The 
cross combinations are indicated at the top of the profiles, respectively. 
G and F indicate parental EG4 and the F, plants, respectively. White and 



black arrowheads indicate the bands representing the de novo mPing 
insertion and the band derived from Nipponbare genome, respectively. 
(B) Mean numbers of de novo mPing insertions in a single F, plant and 
in a self-pollinated plant. The cross combinations are indicated at the 
bottom of the profile. All 16 possible primer combinations were 
analyzed, and mean values were calculated using 1 6 individuals (n = 1 6). 
Bars indicate SE. 

doi:10.1371/journal.pgen.1004396.g001 

Flowering plants have evolved a unique reproductive process 
called double fertilization. In this process, either of two sperm cells 
in pollen fuses with either an egg cell or a central cell in the ovule, 
and then the egg cell fertilized with the sperm cell initiates 
embryogenesis [27]. In rice, the SAM and radicle are regionalized 
in the embryo 3 DAP, and three leaves and the radicle are already 
present in the mature embryo [28]. We detected only three 
radicle/ shoot-specific insertions (Figure 2C), indicating that mPing 
scarcely transposes during the period from the onset of gameto- 
genesis to the early stage (until 3 DAP) of embryogenesis. Among 
the 228 de novo insertions, 1 16 and 17 were shoot-specific and leaf- 
specific insertions, respectively (Figure 2C). This indicates that 
mPing actively transposes in the embryo during the period from the 
regionalization of SAM and radicle (at 3 DAP) to the formation of 
the 3rd leaf primordia (at 8 DAP). Of the 1 33 shoot- and leaf- 
specific insertions, 108 were of the 1st leaf blade (Figure 2D). Since 
the 1 St leaf primordium is formed at 5 DAP, the most active phase 
of the mPing transposition was considered to be from 3 to 5 DAP. 
We detected a large number of radicle-specific insertions as well as 
shoot-specific insertions, and the sum of these insertions accounted 
for 90% of all insertions detected in this study (Figure 2C). Taken 
together, we concluded that mPing in EG4 most actively transposes 
in the 3 to 5 DAP embryo. 

Endosperm is a triploid tissue that is produced by fusing a 
central cell containing two polar nuclei with one of two sperm cells 
in no particular order. The endosperm formation occurs in 
parallel with embryogenesis. The endosperm-specific insertions 
result from the mPing transposition occurred in either gametogen- 
esis or endosperm formation. We observed only two endosperm- 
specific insertions (Figure 2C), supporting that mPing scarcely 
transposes during the period from the onset of gametogenesis to 
the early stage of embryogenesis. The relationship between the 
banding patterns obtained in TD analysis and the timing of mPing 
transposition is summarized in Figure S5. 

Inheritance of de novo mPing insertions to the next 
generation 

In order for mPing to amplify, the de novo insertions must be 
transmitted to the next generation. We performed TD analysis 
using 12 progenies (S2) derived from the main culm and the 
primary tiller of a single selfed parent (Si) to investigate whether 
the de novo insertions detected in oiitogenical analysis are 
inheritable (Figure S6). Both radicle-specific and leaf-specific 
insertions in the Si plants were not detected in the Sj progenies (0 
of 15, 0 of 2, respectively). In contrast, 85% (1 1 of 13) of the shoot- 
specific insertions that were detected in the Si plants also appeared 
in the Sg progenies. This value (85%) is consistent with the 
estimated number of inheritable de novo insertions in our previous 
report [21]. Thus most of the de novo insertions that arose in the 3 
to 5 DAP embryo were successfully inherited to the next 
generation. 

Excisions of mPing during ontogeny of EG4 plants 

We have already determined the sites of all mPing insertions 
(1163 in total) in the EG4 genome [13], and have investigated 
mPing excisions in a small EG4 population using locus-specific 
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Figure 2. De novo mPing insertions during rice ontogeny. (A) Experimental setup for the ontogenical analysis to determine the timing of mPing 
transposition in EG4. Eight progenies (S,) derived from a single parental EG4 plant (So) were grown in a greenhouse. Endosperm, radicle, and leaf 
blades (1 st to 5th) of each S, plant were sampled for DNA extraction. S2 seeds were harvested from the main culm and the primary tiller of each S, 
plant to investigate the inheritance of de novo mPing insertions. The 2nd leaf blade of Sq and S2 plants was also sampled. Detailed information is 
provided in the 'Materials and Methods'. (B) mPing insertions were detected by transposon display. Representative images of shoot-, leaf-, 
endosperm-, radicle-, and radicle/shoot-specific insertions are shown. White arrowheads indicate the bands representing the de novo mPing insertion. 
E: endosperm, R: radicle, L1-L5: 1st to 5th leaf blades. (C) The number of organ-specific de novo insertions in EG4. All 16 possible primer combinations 
were analyzed. En: endosperm-specific insertion, RS: radicle/shoot-specific insertion, R: radicle-specific insertion. Shoot: shoot-specific insertion. Leaf: 
leaf-specific insertion. (D) Percentage of leaf positions where the first de novo mPing insertion was found. Shoot: shoot-specific insertion. Leaf: leaf- 
specific insertion. 

doi:10.1371/journal.pgen.1004396.g002 
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primer pairs [29], [30]. Here we examined the timing of the mPing 
excision with locus-specific PGR using the genomic DNA samples 
that were used for the ontogenical analysis of the de novo insertion. 
We randomly chose 48 markers for this study (Table S3). We 
divided the mPing excisions into five types with the same criteria as 
those used for the de novo insertions: endosperm-, radicle-, leaf-, 
shoot-, and radicle/shoot-specific excisions (Figure S7). There 
were no endosperm-specific and radicle/shoot-specific excisions, 
indicating that no mPing transposition occurs during the period 
from the onset of gametogenesis to the early stage of embryogen- 
esis. We detected seven radicle-specific, six leaf-specific, and three 
shoot-specific excisions. All shoot-specific excisions were detected 
from the 1st leaf blade sample. These results indicate that, lil<e the 
de novo insertion, the mPing excision also occurs during the period 
from the regionalization of the SAM and radicle to the formation 
of the first leaf primordium, although we cannot completely rule 
out the possibility that these excisions occur also in somatic cells of 
mature tissues. Thus, in addition to the experimental results of the 
de novo insertion, we concluded that mPing of EG4 was most actively 
transposing in the 3 to 5 DAP embryo. 

Expression pattern of Ping in EG4 

Both Ping and Pong provide a Myb-like protein and a 
transposase, which are encoded by their ORFl and ORF2, 
respectively (Figure 3A), and have been considered as autonomous 
elements responsible for the mPing transposition. We investigated 
the expression of Hre^-ORFl, Ping-ORF'I, Pong-OKFl, and Pong- 
ORF2 during embryogenesis to evaluate which autonomous 
element plays a predominant role in driving the mPing transpo- 
sition in EG4. Reverse transcription-PGR analysis revealed that 
Ping-OBJ<\ and Ping-GRF2 constitutively expressed in the ovary 
during embryogenesis (Figure 3B). On the other hand, no 
transcriptions of Pong-OKEl and Po«g-ORF2 (Figure 3B) were 
observed. This strongly suggests that Ping predominantly controls 
the mPing transposition in EG4. 

We performed real-time quantitative PGR (qPGR) analysis to 
compare the expression level of Ping-G&Fl and -ORF2 between 
EG4 and Nipponbare during embryogenesis. In all developmental 
stages from 1 to 6 DAP, the expression levels of both ft'n^-ORFI 
and -ORF2 were higher in EG4 than in Nipponbare (Figure 3G, 
D). Since EG4 harbors seven copies Ping, whereas Nipponbare 
has only one copy (Table S4), the difference in the expression 
levels between EG4 and Nipponbare is considered to be 
attributable to the different copy number of Ping. However, we 
found that Ping of EG4 showed different expression patterns from 
that of Nipponbare. In Nipponbare, the expression level of Ping- 
ORFl and -ORF2 gradually dechned until 3 DAP, and restored to 
the basal level at 6 DAP. In contrast, in EG4, the expression levels 
of both Pfng-ORF 1 and -ORF2 rapidly increased, with a peak at 3 
DAP (Figure 3G, D). The ratio of relative expression level (EG4/ 
Nipponbare) clearly demonstrated that Ping might be up-regulated 
in a developmental stage-specific manner in the ovary of EG4 
(Figure 3E). Since mPing transposed during the period from 3 to 5 
DAP, the rapid increase in Ping expression most likely drive the 
mPing transposition. 

Accumulation of Ping transcripts in the embryo triggers 

mPing transposition 

We investigated the spatial expression pattern of Ping by in situ 
hybridization using Ping-speciFic probes. The probe positions were 
indicated in Figure 3A. The Ping transcripts were detected in all 
tissues, viz. embryo, endosperm, and ovary wall, in both EG4 and 
Nipponbare (Figure 4A-G, S8). Among the tissues, the 3 DAP 
embryo of EG4 yielded an exceptionally strong signal, indicating a 



high accumulation of Hng transcripts (Figure 4A), whereas the 5 
DAP embryo showed a much lower accumulation of Hng 
transcripts in EG4 (Figure 4D-F). Such a drastic change in 

accumulation quantity of Ping transcripts with the advance of 
embr\'ogenesis was consistent with the change in the expression 
quantity of Ping with the advance of embryogenesis, which was 
investigated by real-time qPGR (Figure 3G-E). These results 
suggest that the tissue- and developmental stage-specific accumu- 
lation of the Ping transcripts triggers mPing transposition at this 
stage in EG4. To confirm this hypothesis, we evaluated the spatial 
expression pattern of Ping in the SAM during the vegetative phase. 
As described above, mPing hardly transposes in the SAM during 
this phase. The Ping transcripts were detected in all tissues 
including the SAM, and, as expected, there was no obvious 
difference in the signal intensity between EG4 and Nipponbare 
(Figure 4G-I). Thus the Ping transcripts proved to accumulate 
developmental stage-specifically only in the tissue where mPing 
actively transposes. We therefore concluded that the high 
accumulation Ping transcripts triggers the transposition oi mPing 
in the 3 DAP embryo of EG4. 

SNP in an intronic region of P/ng-ORFI 

EG4 has seven Ping elements {Hng-1 to -7), whereas Nipponbare 
has only one (Ping-J\) (Table S4). When we sequenced and 
compared all Ping elements, a single nucleotide polymorphism 
(SNP) in the first intronic region of Ping-OKEl was detected 
between EG4 and Nipponbare (Figure 5A). Ping-M has a 'T' 
nucleotide on the SNP region, whereas all Ping elements in EG4 
have a 'G' nucleotide. We named the former 'T-type f^ng' and the 
latter 'G-type Pin£. 

In addition to EG4, several Aikoku and Gimbozu landraces 
(hereafter AG strains) are known to exhibit high ml^ng activity 
[21]. We investigated the SNP- type of Ping and the copy number 
of Ping and mPing in 93 AG strains, and evaluated the effect of G- 
type Ping on the mPing activity. These 93 AG strains were divided 
into three groups according to the SNP-type of the Ping allele 
(Table S4): strains harboring G-type Ping; strains harboring T-type 
Ping; and strains harboring no Hng. The strains with C-type Ping 
had more mPing copies than those with T-type Ping or no Ping 
(Figure ,5B, Steel-Dwass test, p<0.01). This implies that the C- 
type Ping could drive the mPing transposition. We further 
investigated the expression patterns of /"(Hg-ORFl and -ORF2 in 
two mPing-active strains (A119 and Ai23; and two uiPing-mnctive 
strains (A105 and G190) during embryogenesis (from 1 to 6 DAP). 
A119 and A123 have six and ten copies of G-type Ping, 
respectively, and both A 105 and G190 have one copy of T-type 
Ping (Table S4). Expression analysis revealed that A 105 and G190 
kept low expression levels of ft'ng-ORFl and -ORF2, whereas 
Al 19 and A123 showed high expression levels with a peak around 
3 DAP (Figure 5G-F). This indicates that, in EG4, A119, and 
A123, the developmental stage-specific expression of Ping is 
controlled by the same factor(s) described in the Discussion. 

Discussion 

Our final goal was to elucidate how MITEs attain their high 
copy numbers in the genome. To this end, we chose mPing, which 
is the only active MITE identified in rice, as a material and 
analyzed the timing of mPing transposition in the mPf«g-active 
strain EG4. Gonsequentiy, we successfully found one mechanism 
of the mPing amplification; mPing most actively transposes during 
the period from the regionalization of the SAM and radicle to the 
formation of the first leaf primordium (3 to 5 DAP) by the 
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Figure 3. /'/n^ expression during seed development. (A) Structure of the Ping and Pong elements. Terminal inverted repeats are indicated by 
black triangles. Boxes represent 0RF1 and 0RF2, respectively. Gray horizontal arrows indicate the direction of transcription. RNA probes used are 
indicated below the ORFs. (B) Reverse-transcription PGR analysis of P//ig-ORF 1 , P//ig-0RF2, Pong-ORFI, and Pong-ORF2. Numbers in parentheses are 
PGR cycle numbers. PG: positive control (0.1 ng genomic DNA), NG: negative control (non-reverse-transcribed RNA). (G) Real-time quantitative PGR of 
P/ng-ORFI and (D) P/ng-0RF2. The expression level in the Nipponbare ovary just after pollination was set as 1 . The results are presented as means of 
three biological replicates. Bars indicate SE. (E) The ratio oi Ping expression level of EG4 to that of Nipponbare. The means in (G) and (D) were used for 
calculation. 

doi:1 0.1 371 /journal.pgen.1 004396.g003 



developmental stage-specific up-regulation of the autonomous 
element Ping. 

The transpositions of TEs are categorized into germinal and 
somatic types according to the type of cells where the transposition 
takes place. LOREla in Lotus japonicus is activated in plants 
regenerated from de-differentiated cells and transposes in male 
germ cells by the pollen grain-specific LOREla transcription, 
resulting in the asymmetric transposition of LOREla in the 
reciprocal crosses between the active and non-active lines [31]. 
Tagl in Arabidopsis shows germinal transposition activity in both 
male and female germ cells. Consequently, the reciprocal crosses 



show symmetric transposition of Tagl [32]. These results 
demonstrate that the transposition activity in reciprocal crosses 
reflects the tissue specificity of germinal transposition. In this study, 
reciprocal crosses between EG4 and Nipponbare showed the same 
mPing transposition activity, which may suggest that mPing in EG4 
transposes in both male and female germ cells. However, we 
obtained only a few de novo endosperm-specific and radicle/ shoot- 
specific insertions in the ontogenical analysis, although we detected 
a number of de novo shoot-specific and radicle-specific insertions. 
We therefore concluded that most mPing transposes not in germ 
cells but in somatic cells after pollination. Somatic transposition 
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Figure 4. Detection of Ping-ORF^ spatial expression patterns by in situ liybridization analysis. Longitudinal sections through the ovary 3 
days after pollination of (A, C) EG4 and (B) Nipponbare; the ovary 5 days after pollination of (D, F) EG4 and (E) Nipponbare; and the shoot apical 
meristem of (G, I) EG4 and (H) Nipponbare seedlings were hybridized with antisense (A, B, D, E, G, H) or sense (C, F, I) RNA probes. Little staining was 
obtained with the sense probe (F). Em: embryo, En: endosperm, OW: ovary wall, SAM: shoot apical meristem. Scale bars represent 100 |im. 
doi:10.1371/journal.pgen.1004396.g004 



that occurs at the late stage of plant development often produces 
spotted and striped segments in tissues, such as maize seed coat 
variegation caused by Matotor excision from the bz2 gene [33], [34] 
and rice leaf color variegation by aDari/ excision from the OsClpPS 
gene [35] . In animals, somatic transposition is seldom transmitted 
to the next generation because germ cells are set aside from 
somatic cells at the early stage of embryogenesis. On the other 
hand, in plants, germ cells are generated from somatic cells at the 
reproductive stage. In rice, gametes are generated in the SAM; 
therefore, somatic transposition that occurred in the SAM can be 
transmitted to the next generation via gametes. In this study, we 
revealed that most mPing elements transposed in somatic cells of 
the embryo during the period from 3 to 5 DAP. Being a class II 
TE that transposes by a cut-and-paste mechanism, mPing is 
expected to be eliminated from genomic DNA with a certain 
frequency. However, a previous report demonstrated that the 
mPing excision sites would be repaired by utilizing a copy of mPing 
from either the sister chromatid or from the homologous 
chromosome [29]. The mPing excision site cannot be repaired if 
mPing transposes in germ cells, which are haploid. It is therefore 
considered that the somatic transposition of mPing is an important 
factor for mPing amplification in the genome. 

The autonomous elements Ping and Pong mediate mPing 
transposition in the rice genome, yiaay japonica cultivars, including 
EG4 and Nipponbare, have both Ping and Pong. This study 
demonstrated that Ping plays a predominant role in mPing 
transposition in EG4. However, a heterologous expression assay 
using Arabidopsis and yeast showed that Pong had a higher 
catalytic capacity for mPing transposition than Ping [22], [23]. 
Furthermore, mPing transposition was observed under stress 
conditions in several rice cultivars harboring only Pong [14], 
[17], [19]. In this study, however, we detected very low expression 



of Pong through the development of rice plants, indicating that Pong 
would be epigeneticaUy silenced at the transcriptional level in 
EG4. In contrast. Ping constitutively expressed in aU organs 
including the SAM and embryo. Nevertheless, mPing could be 
transposing most actively in the embryo during the period from 3 
to 5 DAP. Since the stage-specific up-regulation of Ping was 
observed during the period of mPing transposition, we hypothe- 
sized that the expression level of Ping needed to exceed a certain 
threshold of mPing transposition. 

AU mPm^-active strains (EG4, Al 19, and A123) showed higher 
expression of Ping with a peak around 3 DAP than the mPing- 
iiiactive strains (Nipponbare, A 105, and G190). Although further 
experiments are needed to elucidate the mechanism of develop- 
mental stage-specific up-regulation of Ping expression, we propose 
two hypotheses: (1) position- and dosage-effect, and (2) effect of 
SNP. The details of the hypotheses are as follows. 

Position- and dosage-effect 

Chromosomal position and copy number of TE often affect the 
transposition activity. The former is known as 'position effect' and 
the latter as 'dosage effect'. Eight independent Tam3 copies 
residing in the Antirrhinum majus genome show different transpo- 
sition activities from each other [36]. In Arabidopsis, germinal 
reversion frequency of Tagl increases in proportion to its copy 
number [32]. The mPm^-inactive strains Nipponbare, A 105, and 
G190 have only one Ping at the same locus, whereas the mPing- 
active strains EG4, A119, and A123 have respectively seven, six, 
and ten copies of Ping at different loci except for the Ping-1 locus. 
Furthermore, the expression pattern oiPing showed slight variation 
among the mPm^-active strains harboring only C-type Ping. These 
results suggest that the developmental stage-specific up-regulation 
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Figure 5. SNP in tlie first intronic region of Ping-ORFl. (A) Determination of the SNP sequence in the first intronic region of P/ng-ORFI. The 
arrowhead indicates the position of the SNP. The number indicates the position of the Ping element. Ping harboring +1261C SNP and +1261T are 
named 'C-type' and 'T-type' Ping, respectively. (B) Box plots of mPing copy number in AG lines. The top and bottom of the boxes mark the first and 
third quartiles, respectively. The center line represents the median, and the whiskers show the range of observed values within 1.5 times the 
interquartile range from the hinges. Values beyond 1 .5 times the interquartile range from the nearest hinge are marked by open circles. 'No Ping', 'C- 
type Ping,' and 'T-type Ping' indicate the groups having no Ping, C-type Ping, and T-type Ping, respectively. Expression of (C) P/ng-ORFI and (D) Ping- 
0RF2 during embryogenesis in mP/ng-active strains (A1 19 and A123) and mP//ig-inactive strains (A105 and G190). The results are presented as means 
of three biological replicates. Bars indicate SE. The ratio of (E) P/ng-ORFI and (F) -0RF2 expression level of A105, A119, A123, and G190 to that of 
Nipponbare. The means in (Fig. 3C) and (Fig. 3D) were used for calculation. 
doi:10.1371/journal.pgen.1004396.g005 
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oiPing expression is probably regulated by the position-effect and/ 
or the dosage-effect. 

Effect of SNP 

Intronic SNPs are known to cause drastic effects on gene 
expression. In humans, an intronic SNP in SLC22A4 affects 
transcriptional efliciency in vitro, owing to an allelic difference in 
affinity to the transcriptional factor RUNX 1 [3 7] . Furthermore, a 
SNP located in the intronic enhancer region of the thyroid 
hormone receptor P gene enhances pituitary cell-specific tran- 
scriptional activity [38] . In this study, we demonstrated that a SNP 
is present in the intronic region of Pj'ng-ORFl, and Ping elements 
in the AG strains were categorized into either T-type or C-type 
Ping according to the SNP-t}'pe. Since all strains that showed a 
peak in the expression analysis had only C-type Ping, the intronic 
SNP might influence the developmental stage-specific up-regula- 
tion of Pz'reg expression. T-type Ping was. present in 14 AG strains as 
one copy, and its chromosomal location did not differ among 
strains. In contrast, the copy number of C-t>pe Ping varied from 
one to ten, and their chromosomal locations, except for Ping-l, 
differed from each other. These results indicate that T-type Ping 
has lost its activity, whereas C-type Ping may be still active in the 
rice genome. Furthermore, we found that the copy number of 
mPing was significantiy larger in strains harboring C-type Hng than 
in strains harboring T-type Hng. This strongly supports that C- 
t)pc SNPs in the intronic region of Ping contribute to the 
amplification of mPing, presumably by the developmental stage- 
specific up-regulation of Ping expression. 

Since the transposition of TEs often damages the host genome, 
TEs with high transposition activity are targeted by the silencing 
mechanisms. Nevertheless, MITEs amplify to very high copy 
numbers not only in plant genomes but also in animal genomes. 
Very littk' is launMi about h()^\- MITEs attain their high copy 
numbers by escaping the silencing mechanism. The transposition 
of mPing is transiently induced by various stresses [14]-[18], 
indicating that the activity of mHng is suppressed by the silencing 
mechanisms in many cultivars. Thus, mPing must overcome the 
silencing mechanism in order to maintain the transposition activity 
under natural growth conditions. Our results revealed that mPing 
in EG4 ^\-as moljilizcd by the sufficient supply oi Ping transcripts 
produced only during the period otmPing transposition. This stage- 
specific activation is thought to be a strategy of the mPing family to 
amplify mHng by escaping from the silencing mechanism of the 
host genome. Since no active MITEs other than mPing so far have 
been identified, it is very difficult to elucidate if the other MITEs 
also attain their high copy numbers in the same way as mPing 
amplifies. Given that the other active MITEs are identified, 
however, our study will help to understand their amplification 
mechanisms. Our previous study documented the generation of 
new regulatory networks by a subset of mPing insertions that render 
adjacent genes stress inducible [13]. In addition to mPing, other 
MITEs also contribute to gene and genome evolution via 
providing new promoter regulatory sequences, transcriptional 
termination elements, and new alternative exons [39], suggesting 
that the amplification of MITEs causes gene and genome 
evolution. Our results provide clues to further understand not 
only the amplification mechanism of MITEs but also the co- 
evolution of MITEs and the host genome. 

Materials and Methods 

Plant materials and sampling 

EG4 (cultivar Gimbozu), Nipponbare, and 94 Aikoku/ Gimbozu 
landraces were used in this study (Table S4). Aikoku/Gimbozu 



landraces were provided from the GenBank project of the 
National Institute of Agrobiological Science, Ibaraki, Japan. 
Reciprocal crosses between EG4 and Nipponbare were made in 
a green house. Before pollination, all anthers were removed from 
the flowers of maternal plants. The pollinated flowers were 
covered with protective bags to prevent outcrossing until harvest. 
After harvesting, success of crosses was checked with the molecular 
markers. For ontogenical analysis, eight progenies of EG4 (Si) 
derived from a single parental plant (Sq) were grown in a 
greenhouse, and all S2 seeds were harvested. For S] plants, each 
seed was cut into two halves, and the half including the embryo 
was germinated and the other was sampled. After germination, the 
radicle and the 1st, 2nd, 3rd, 4th, and 5th leaf blades were 
sampled. The second leaf was collected from Sq and S2 plants. For 
estimation of Hng and mHng copy numbers, eight bulked plants 
were sampled. For RNA extraction, ovaries before pollination and 
ovaries at 1,2, 3, 4, 5, and 6 DAP were collected. AU samples were 
immediately frozen in hquid nitrogen and stored at — 80°C until 
use. 

DNA extraction and transposon display 

DNA extraction and transposon display was performed 
following a published protocol [30]. For DNA extraction from 
endosperm, we used GM quicker 2 (Nippon Gene). 

Locus-specific PCR 

Sequencing of mPfn^-flanking fragments excised from transpo- 
son display gels and primer design were performed following a 
published protocol [30]. The genomic locations of the mPing 
insertion sites were forecasted by a BLAST search in the Rice 
Annotation Project Database (RAP-DB; http://rapdb.dna.affrc. 
gojp/) [40], [41] using mPing flanking sequences as queries. To 
prepare enough templates for PCR, whole genome amplification 
was performed using an Ulustra GenomiPhi V2 Kit (GE 
Healthcare). mPing excision was detected by PCR with mPing- 
sequence characterized amplified region (SCAR) markers [29]. 
PCR was performed in 10-|il reaction volumes containing 10 ng of 
the template DNA, 5 |il of GoTaq Green Master mix (Promega), 
5% DMSO, and 0.25 nM of each primer (Table S3). PCR 
conditions were as follows: 94°C for 3 min; 40 cycles of 98°C for 
10 s, 57°C for 30 s, and 72°C for 45 s; and 72°C for 5 min. To 
detect the presence oi Ping-N, -1, -2, -3, -4, -5, -6, and -7, eight 
Pfrag-SCAR markers were used. The genomic locations of the Ping 
insertion sites were referred from a previous report [42]. For 
detection of the Ping-l allele, PCR was performed in 10-|J.l reaction 
volumes containing 10 ng of template DNA, 0.2 U of KOD FX 
Neo (Toyobo), 1 xPCR buffer for KOD FX Neo (Toyobo), and 
0.2 |xM of each primer (Table S5). PCR conditions were as 
follows: 94°C for 3 min; 35 cycles of 98°C for 10 s, 60°C for 30 s, 
and 68°C for 90 s; and 72°C for 1 min. For detection oiHng-N, -2, 
-3, -4, -5, -6, and -7 alleles, PCR was performed in lO-^ll reaction 
volumes containing 1 0 ng of template DNA, 5 |il of GoTaq Green 
Master mix (Promega), 5% DMSO, and 0.25 |xM of each primer 
(Table S5). PCR conditions were as follows: 94°C for 3 min; 35 
cycles of 98°C for 10 s, 60°C for 30 s, and 72°C for 45 s; and 
72°C for 1 min. 

RNA isolation and expression analysis 

Total RNA was isolated using TriPure isolation reagent (Roche) 
and digested using RNase-free DNase (TaKaRa). First strand 
cDNA was synthesized using a Transcriptor first strand cDNA 
synthesis kit (Roche). For reverse transcription PCR, PCR was 
performed in 10 |ll reaction volumes containing cDNA generated 
from 4 ng total RNA, 0.2 U of KOD FX Neo (Toyobo), 1 xPCR 
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bufier for KOD FX Neo (Toyobo), and 0.5 nM of each primer. 
PGR conditions were as follows: 94°C for 3 min; 35 or 45 cycles of 
98°C for 10 s, 60°C for 10 s, and 68°C for 10 s. Relative 
quantification of Ping-ORF 1 and Ping-OBJ^2 were calculated by 
the 2 ^^'-^^ method [43] using Light cycler 1.5 (Roche). The 
UBQ5 gene was used as the calibrator gene. The thermal profile 
consisted of 10 min at 95°C; and 45 cycles of 4 s at 95°C, 10 s at 
60°C, and 1 s at 72°C. Amplification data were collected at the 
end of each extension step. The primer pairs used in this study are 
listed in Table S6. 

Paraffin sectioning and in situ hybridization 

Plant samples were fixed with 4% (w/v) paraformaldehyde and 
1% Triton X in 0.1 M sodium phosphate buffer for 48 h at 4°C. 
They were then dehydrated in a grad("d c'tlianol series, substituted 
with 1-butanol, and embedded in Paraplast Plus. The samples 
were sectioned at 8-|im thickness using a rotary microtome. 
Fragments of ft'wg-ORFl (1091 bp) and Ping-OUF2 (1368 bp) 
were cloned into pBlueScript SK+ (Stratagene) and sequenced. 
For digoxigenin-labeled antisense/ sense RNA probe .synthesis, in 
vitro transcription was performed using T7 RNA polymerase and 
T3 RNA polymerase. In situ hybridization and immunological 
detection with alkaline phosphatase were performed according to 
Kouchi and Hata [44]. 

SNP detection 

PGR was performed in 10-|J,1 reaction volumes containing 10 ng 
of template of DNA, 5 |il of GoTaq Green Master mix (Promega), 
5% DMSO, and 0.25 nM of each primer. PGR conditions were as 
follows: 94°G for 3 min; 35 cycles of 98°G for 10 s, 60°C for 30 s, 
and 72°G for 30 s; and 72°G for 1 min. PGR primers used in this 
study are listed in Table S6. Because the original sequence 
contained an Afa I restriction site, one mutation was introduced 
into the reverse primer. The 5-^1 PGR products were mixed with 
5 nl restriction mixture containing 1 U Afa I (TaKaRa), 33 mM 
Tris-acetate (pH 7.9), 10 mM Mg-acetate, 0.5 mM Dithiothreitol, 
66 mM K-acetate, and 0.01% (w/v) bovine serum albumin. After 
16 h incubation at 37°G, DNA gel electrophoresis was performed. 
PGR products (502 bp) including +126 IT SNP were not digested, 
whereas PGR products including +1261G SNP were digested into 
two fragments (352 bp and 150 bp). 

Estimation of Ping and mPing copy number 

To determine the copy number ofPinghy Southern blot analysis, 
genomic DNA samples were digested with Eco RI restriction 
enzyme. These samples were loaded onto an agarose gel, separated 
by electrophoresis, blotted onto a nylon membrane, and probed 
with the Ping fragment. The mPing copy number was determined by 
real-time quantitative PGR as described previously [45] with litde 
modification. Quantitative PGR was performed using the Light- 
Gycler 480 system (Roche). PGR was performed in 20 |tl reaction 
volumes containing 5 Jil genomic DNA (0.4 ng/^ll), 1 xlightCycler 
480 SYBR Green I Master mix (Roche), and 0.5 jlM of each 
primer. Specificity of the amplified PGR product was assessed by 
performing a melting curve analysis on the LightGycler 480 system. 

Supporting Information 

Figure SI Verification of the reciprocal crosses between EG4 
and Nipponbare. (A) Locus-specific PGR analysis of Ping with 
locus-specific markers. The maker names and cross combinations 
are indicated at the top of the profiles. (B) Locus-specific PGR 
analysis of mPing in Nipponbare genome. The genomic location of 
the mPing insertions and cross combinations are indicated at the 



top of the profiles. G and F indicate parental EG4 plants and the 
Fi plants, respectively. Lane M: DNA size marker (Gene Ladder 
100, Nippon Gene), Lane N: Nipponbare. Black and white 
arrowheads show the bands indicating the presence and absence of 
Ping/mPing, respectively. 
(TIF) 

Figure S2 Locus-specific PGR analysis oi de novo mPing insertions 
in Fi progenies. Ten representati\e results are shown. The 
genomic locations of the mPing insertions and cross combinations 
are indicated at the top of the profiles. G and F indicate parental 
EG4 and the Fj plants, respectively. Lane M: DNA size marker 
(Gene Ladder 100, Nippon Gene), Lane N: Nipponbare. Black 
and white arrowheads show the bands indicating the presence and 
absence of mPing, respectively. 
(TIF) 

Figure S3 Ontogenical analysis of mPing transposition in EG4 by 
transposon display. Eight progenies (Sj) were derived from a single 
parental EG4 plant (So). The 2nd leaf blade of the So plant and the 
endosperm, radicle, and leaf blades of each Si plant were sampled 
and subjected to transposon display. White, red, and green 
arrowheads indicate shoot-, radicle-, and leaf-specific insertions, 
respectively. The rice plant has alternate distichous leaves; 
therefore, we analyzed the insertion in both [n+l]th and [n+2]th 
leaves to confirm whether the insertion detected in [n]th leaf is 
leaf-specific or shoot-specific. But we did not investigate the 
specificity of the insertions detected in the 4th and 5th leaves using 
their upper leaves. For this reason, we did not categorize such 
insertions and marked with the gray arrowhead. E: endosperm; R: 
radicle; L1-L5: 1st to 5th leaf For progeny 2-8, samples are 
applied in the same order as for progeny 1 . 
(TIF) 

Figure S4 Locus-specific PGR analysis of de novo mHng insertions 
in various tissues of a single EG4 plant. Eight representative results 
are shown. The genomic locations of the mPing insertions and 

insertion types are indicated at the top of the profiles. Lane M: 
DNA size marker (Gene Ladder 100, Nippon Gene), Lane N: 
Nipponbare, Lane E: endosperm. Lane R: radicle, Lane L1-L5: 
1st to 5th leaf Black and white arrowheads show the bands 
indicating the presence and absence of mFHng, respectively. 
(TIF) 

Figure S5 Schematic representation of the relationship between 
banding patterns obtained in transposon display and the timing of 

mPing transposition. If mPing transposes in the period indicated by 
the red bar, the schematic banding patterns indicated by the 
arrows will be obtained. E: endosperm, R: radicle, L1-L5: 1st to 
5th leaf blade, DAP: days after pollination, SAM: shoot apical 
meristem. 
(TIF) 

Figure S6 Inheritance of de novo mPing insertions in EG4. S2 
plants derived from the main culm and the primary tiller of a 
single Si plant were assayed. The shoot-specific insertion in the Si 
plant (white arrowhead) was inherited by S2 plants, whereas the 
radicle-specific insertion (red arrowhead) was not. E: endosperm; 
R: radicle; L1-L5: 1st to 5th leaf 
(TIF) 

Figure S7 mPing excisions in EG4. mPing excisions were detected 
by locus-specific PGR using the genomic DNA samples that were 
used for the ontogenical analysis of the de novo insertion. We 
analyzed 48 loci. Black and white arrowheads show the bands 
indicating the presence and absence of mPing, respectively. Figures 
indicate (A, B) shoot-specific excisions, (G, D) leaf-specific 
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excisions, and (E) radicle-specific excision. G: EG4 (Sq plant); N: 
Nipponbare; E: endosperm; R: radicle; L1-L5: 1st to 5th leaf. 
(TIF) 

Figure S8 Detection oiPing-OK¥2 spatial expression patterns by 
in situ hybridization analysis. Longitudinal sections through the 
ovary 3 days after pollination of (A) EG4 and (B, C) Nipponbare 
were hybridized with (A, B) antisense or (C) sense RNA 
probes. 
(TIF) 

Table SI 

(XLSX) 



De now mPing insertion sites detected in Fi progenies. 



Table S2 De novo mPing insertion sites detected in various tissues 
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